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Q . ABSTRACT 
<D 

f"*^ ■ We propose simple, explicit, analytical approximations for the kinematics of an adiabatic 

blast wave propagating in an exponentially stratified ambient medium, and for the onset of 
radiative cooling, which ends the adiabatic era. Our method, based on the Kompaneets implicit 
solution and the Kahn approximation for the radiative cooling coefficient, gives straightforward 
estimates for the size, expansion velocity, and progression of cooling times over the surface, when 
applied to supernova remnants (SNRs). The remnant shape is remarkably close to spherical for 
moderate density gradients, but even a small gradient in ambient density causes the cooling time 
to vary substantially over the remnant's surface, so that for a considerable period there will be a 
cold dense expanding shell covering only a part of the remnant. Our approximation provides an 
effective tool for identifying the approximate parameters when planning 2-dimensional numerical 
models of SNRs, the example of W44 being given in a subsequent paper. 
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Introduction 



x 

A good model of a supernova remnant's (SNR's) blast wave propagating through the ambient medium 
allows us to derive, from the observed quantities, important properties of the medium, as well as the 
explosion energy. The Sedov solution (1959) for adiabatic expansion into an isotropic medium usually serves 
as the prototype, and further, more realistic models also assume isotropy (see e.g. Jun, Jones & Norman 
1996). Nevertheless, in most cases of astrophysical interest (stellar winds, Garcia-Segura & Mac Low 1995; 
SNRs in vicinity of molecular clouds, Dohm-Palmer & Jones 1996; galactic supershells, Maciejewski et al. 
1996; and starburst galaxies, Jogee, Kenney & Smith 1998), the ambient medium is not isotropic. At the 
next level of approximation, Kompaneets (1960) found an implicit solution for a blast wave propagating in 
an exponentially stratified ambient medium, under the assumption that the post-shock pressure is uniform 
over the surface of the remnant. An explicit form can be obtained when one additional assumption is 
made, that the shape of the shock wave can be approximated by a prolate ellipsoid. After an alternative 
presentation of the Kompaneets solution in §2, we derive explicit analytical formulae for the size, velocity 
and the onset of cooling in §3. A prolate ellipsoid is a good approximation for all evolutionary times prior to 
blowout (Kompaneets 1960), but we show in §4 that the Kompaneets approximation of uniform post-shock 
pressure loses validity much earlier. In §5 we show how our explicit approximation can be used in exploring 
the parameter space to identify values of interest in making 2-dimensional numerical models of SNRs. 
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2. A simple method to derive basic Kompaneets conclusions 

Before we introduce our assumption, let us show a simple way of re-deriving the basics of the original 
Kompaneets solution, using one of his conclusions, namely that in the direction perpendicular to the density 
gradient, the remnant is widest halfway between the high and low ambient density extremes of the shock, 
a location we shall refer to as the equator. With the resulting formulae for the remnant size parallel and 
perpendicular to the gradient, our approximation for the shock shape immediately provides the volume and 
pressure needed to learn the time evolution. 

We consider the shock expansion in cylindrical coordinates (z, r): the explosion occurs at (z, r) = (0, 0) 
(Fig.l). The ambient density po is an exponential function of z: 

Po = P*e- Z '\ (1) 

where is the density at the explosion site and h is the stratification scaleheight. Denoting zh and zj, as 
the location of the dense and the tenuous end of the remnant, respectively, we have z^ — zh = 2a, where 2a 
is the extent of the remnant in the direction along the density gradient. For a strong, non-radiative shock 
in a 7 = 5/3 gas, the post-shock pressure Ps, pre-shock density p n and shock velocity v are related by 

Pov 2 = ^Ps- (2) 

Assuming after Kompaneets that the post-shock pressure Ps is constant over the surface of the remnant, 
the shock speed at a particular time can be written as a function of z : 

v(z,t)=v*(t)e z / 2h . (3) 

Comparing the values of shock velocities at the ambient density extremes and at the equator, we can 
eliminate the unknown v* (t) and calculate the evolution of the remnant shape as a function of its semi- major 
axis a. 

The rate of expansion along the density gradient is an average of shock velocities at two density 
extremes: vr and vl 

. V H +V L e z H /2h + e z L /2h 

a= g = V *^ 2 ' () 

where vr and vl are velocity magnitudes, zh is negative, and vh = —zh- The center of the remnant is 
displaced from the explosion site by s — 0.5(zl + zh)- At z = s (the equator), the remnant expands only in 
the lateral direction and the expansion rate is 

v s =v.(t)e'' 2h . (5) 

Although the equator moves along the z axis (Fig.l) as the remnant expands, the rate of change of the 
remnant's semi-minor axis b is exactly equal to v$: b = v$- The evolution of the displacement s can be 
written as 

. ZL+ZH V L -V H e z L /2h _ e z H /2h 

*=^— = =«•(*) ~ 2 • ( 6 ) 

After substituting zr = s — a and zl = s + a and simple manipulations, we can eliminate the constant v*(i) 
and get 

da a ds a 

— = cosh— ; — =tann— . (7) 

db 2h ' da 2h w 



- 3 - 




Fig. 1. — (left) Representation of SNR characteristics used in the text. The ellipsoid is axially symmetric 
about the z axis. Fig. 2. — (right) The ratio of cooling times on the tenuous and dense ends of the remnant 
expanding in a stratified medium, as a function of the ratio of extreme ambient densities of an equivalent 
remnant at its cooling time. The equivalent remnant is defined as a remnant with the same explosion energy, 
but expanding in a uniform density equal to the density at the explosion site. For the radius r coo i of the 
equivalent remnant at its cooling time, the ratio of extreme ambient densities is defined as W — exp(2r coo \/h), 
where h is the density scaleheight of the ambient medium. 



Integrating the equations above, we can readily find b and s as functions of a 

b . , a s , a /0 . 

tan— = sinn — ; exp — = cosh — -. (8) 

2h 2h * 2h 2h w 

These results are independent of the assumed shape of the blast wave and the same as found by Kompaneets 
(1960), after his auxiliary variable x is recognized as tanh(a/2/i). In particular, the equations above allow 
us to estimate the flattening of the remnant b/a and the relative displacement of the center s/a. It appears 
that even remnants extending several scale-heights retain a nearly spherical shape: for example, if a = 1.7h, 
the ambient density contrast between the ends of the shock is almost 30, but the remnant remains barely 
flattened (size ratio 6/a=0.9). On the other hand, the ellipsoid's center is displaced from the explosion site 
by s = 0.38a in the low density direction. The same conclusions about sphericity and the center shift can 
be drawn from the analytical solution based on sectoral approximation (Gnatyk 1988). Hydrodynamical 
models constructed by Dohm-Palmer & Jones (1996) arrive to these conclusions as well. One can also 
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notice that exp(zf//2/i) = exp((s — a)/2h) — 0.5 (1 + exp(— a/h)): in the blowout case, when a — > +oo, 
it is 0.5. This is Kompaneets' famous result, that as blowout occurs and Zl becomes large, zh goes to a 
constant —2ft, log 2. The lateral size at blowout also agrees with Kompaneets' prediction: for a — > +oo, we 
get b/2h -> tt/2. 



3. Our approximations 

The implicit solution by Kompaneets shows that the shock has the shape of a somewhat boxy ellipsoid. 
In fact, its shape is virtually indistinguishable from a true ellipsoid within limits of accuracy set by 
Kompaneets assumption of uniform post-shock pressure. We use this fact in order to introduce time as an 
independent variable instead of a, when we assume that the post-shock pressure Ps has the same volume 
relationship as in the ID Sedov solution (see e.g. Bisnovatyi-Kogan & Silich 1995) 

P s = ^~ (9) 
3 25 V 1 ' 

where £ 5 = 2.025, E is the explosion energy and V = ^irab 2 is the volume of the remnant approximated 
by an ellipsoid. Combining equations (|I|), (||) and ||), and identifying z — s along the minor axis of the 
ellipsoid we have 

p„ e -.M \? = JL. (10) 



After expressing s and b as functions of a, we get 



which, after integrating gives 




25 ab 2 ' 



cosh 2 (a/2/i) 



a arctansinh(a/2ft) ' 



E 
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t(a) = C 5/2 f — ) (2ft) 5/2 I(a/2h) = t 2h I(a/2h), (12) 

where tih is the nominal time at which a Sedov remnant in a homogeneous medium would reach a radius of 
2h, and the integral I(x) is defined in Table 1 for numerical evaluation. The numerical form of t^h is 

t 2h = 17.32 yr (2/i pc ) 5 / 2 nj* / (13) 

where h pc is h in parsecs, £51 is E in units of 10 51 ergs, and n* is the nuclear number density of the ambient 
medium in cm~ 3 . In this notation, the mass density is p = ran, where the average mass per nucleus is 
m = (1.4/1. l)mjy, run being the mass of hydrogen. 

The shock velocity at any z can be easily calculated by combining equations (||) and (||) and 
substituting the local ambient density for po- I n the particular case of top, bottom and equator 
expansion velocities, one can use equation (^J), which after combining with equation (^) gives 
a = Vjj exp(a/2h) cosh(a/2/i) = vl exp(— a/2h) cosh(a/2h). After substituting a from equation (pd|), 
one gets explicit formulae for shock velocities at two density extremes: vh and v^. Obviously, the shock 
expansion velocity on the equator is the geometrical average of vh and Uj, — these results are given in 
Table 1. The ratio of expansion velocities at the two ends is exp(a/h). 

In SNR evolution in a uniform medium, the adiabatic era is brought to close by the onset of significant 
radiative cooling, at remnant radius r coo \(E, /?*), closely followed by formation of a dense shell. When a 
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remnant is evolving in a density gradient of scaleheight h, cooling occurs first on the dense end, and then 
spreads over the remnant surface. The evolution has two characteristic scales: the kinematic h and the 
radiative r coo i, with a history of shell formation that is dependent on their ratio. 

In order to estimate cooling times, we used Kahn's cooling law approximation (1976), in which the 
cooling coefficient takes the form L = aT^ 1 / 2 . It is a reasonable approximation to the actual cooling 
coefficient for the temperature range of cooling SNRs (Smith et al.1996), and it has the special property 
that in the absence of thermal conduction, the time for a parcel of hot gas to cool from an initial state is 
independent of its history and equals 

At coo] = -gL (14) 

For a newly shocked parcel of gas, the post-shock values of pressure Ps, temperature T, and density n are 
known, so we can find At coo i- This incremental cooling time Ai coo i is added to the time at which a parcel 
is shocked, to find the actual cooling time of that parcel, and then the sum minimized over all parcels to 

3/2 

find the earliest parcel to cool. For a homogeneous medium, Ai coo i is proportional to P s and therefore to 
r -9/2 or £-9/5 Numerically, the minimum of t + Ct~ 9 / 5 occurs at t\ = [9C/5] 5 / 14 such that the first gas to 
cool was initially a distance 

ri = r(t x ) = 19.62 pc E 2 5 { 7 /n 3 J 7 (15) 

from the explosion site, and was shocked at time t\ (see Cox & Anderson 1982, Cox 1986). The expression 
for C can be obtained by substituting equation (||) to (|T3), and using formulae for the Sedov radius and 
velocity from Table 1. One can get then C = ± (Cv/3/5) [m 9 / 10 /(a^/xk)] (^ a6 /nJ- 6 ) . At the minimum, 
Atcooi = for other radii, Ai coo i = |ii(ri/r) 9 / 2 . First parcel's cooling is complete at i coo i = ^ti, when 
the remnant radius is r coo i = (if) 2 ^ 5 ^i. We define t coo \ as the cooling time of the remnant. 

In our approximation for the remnant in a density gradient, a combination of equations ([!]), (j^) and 
(13) gives 

, , 5 , , / rf \ 3/2 fnA 5/2 in El{ 2 exp(5z/2h) , N 

where we substituted a = 1.3 x 1CP 19 K 1 / 2 cm 3 erg s _1 in Kahn's formula for cooling. Noting that 
t 2h = i 1 (2/ l /ri) 5/2 yields 

5 f r i\ 7 exp(5z/2h) , . 

At cooi (a,z) = -t 2h [-) {a/2hmb/2h)3 - (17) 

A measure of the likehood of finding a remnant with a cold shell only on the dense end is provided by the 
ratio t~ ool /t^ ool of cooling times between the tenuous and dense ends. These cooling times can be obtained 
by minimizing the value of the sum At coo i + t(a) of times given by equations (17) and ( |l2| ) over x = a/2h 



tfool = *2h min^ 



(18) 



where the function g±(x) is given in Table 1, and the upper (lower) sign corresponds to the dense (tenuous) 
end. The formula to derive values of x at minimum, xf , is given in Table 1. One can see that x^r is a 
function of only one variable: n/2h = (14/9)~ 2 / 5 r coo i/2ft,. 

The ratio of cooling times between the tenuous and dense ends of the remnant is numerically close 
to the density contrast between the two ends, W — exp(2r coo i//i), at the nominal cooling radius. Figure 
2 displays the ratio of cooling times as a function of W. Note that for cooling time ratios less than 10 
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(essentially all for which the parameter is interesting), it is equal to the density contrast W within 20%. 
From Figure 2 we see that cooling takes more than twice as long on the tenuous end for W > 2.4, or roughly 
speaking, for r coo i > h/2, yielding the unsurprising result that when the density differential between the 
two ends is large at the nominal onset of cooling, the evolution will contain a long period with a partial 
shell. This result may contribute to the explanation of the fact that we usually can see only one side of the 
expanding HI shells (Heiles 1979, Koo & Heiles 1991). 

For completeness, we now provide an approximate evaluation of the time at which the first parcel 
on the equator cools, as a measure of the proximity of a remnant to having half a complete shell. This 
is somewhat more difficult, because the location of the equator is shifting with z. Any parcel cooling 
at the equator was not on the equator when shocked. Nevertheless, we are able to provide a plausible 
approximation for the nearly spherical case. We use equation (^|), which relates the post-shock pressure to 
the explosion energy and volume, but approximate the volume by V = |7ra 3 , and then substitute for a the 
solution of the Sedov ID problem a = (E / p*) 1 / 5 ^ 2 / 5 . The approximate post-shock pressure is then 

Ps = ^C 2 (^ 2 ^) 1/5 *- 6/5 - (19) 
25 

From equations (|lj) and (|l^), the cooling time for a parcel on remnant's minor axis, which is hit by the 
shock at the time t s , when the center of the remnant was at a distance s from the explosion site, is 

EJ3/5 

At coo i(s) cx -— tj 9 ' 5 exp(5s/2/i). (20) 

p* 

By assuming that At coo i(s) rs |t s , which is an exact equality for the uniform medium only, an approximate 
formula for t s can be derived. The corresponding approximate expression for t coo i is 

14 E 3 ^ 14 25 s 

t coo i = Atcooi(s) + t a = —t s = 4.60 x 10V -^j- exp( — -). (21) 

and approaches the cooling time for Sedov solution when h — ► +oo. The dependence on s can be removed 
by expressing s in terms of a from equation (^|) and substituting a from the Sedov ID solution. Thus we 
find the iterative solution for the cooling time on the minor axis 

p 3/14 / /F v 1/5 ( 9_. N2/5' 

■)4, rT . ^51 r^cli 25 / 14 I met 51 I ^ 14 c ° o1 ^ 



tcooi = 4.60 x 10 4 yr cosh^ 14 0.16 V14 ' ) . (22) 



4. Accuracy of the Kompaneets Approximation 

The assumptions of uniform post-shock pressure and an exponential atmosphere lead directly to an 
equation (||) for the displacement s of the remnant center from the explosion site, as a function of the 
semimajor axis, a. In the hydrodynamical models performed for the particular case of SNR W44 by Shelton 
et al.(1998, Paper II), we noticed that this equation appeared to overestimate the displacement of the 
center. That implies that the post-shock pressure at the high density end is higher than at the low density 
end. In the following, we provide a rough justification for such a differential. Note that numerical studies of 
s(a) may yield a reasonably useful test case for intercomparison of 2D hydrocodes. 

The upturn in pressure at the edge of the Sedov remnant can be attributed to the effective gravity 
in the decelerating post-shock gas. The post-shock mass velocity is V2 = jv, so the outwardly directed 
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efFective gravity is roughly — ~v. Assuming that the post-shock gas is in hydrostatic equilibrium, we get 
^ ~ — pg ~ — 3po^j given p = 4po (f? is the coordinate locally perpendicular to the shock front and directed 
outwards). The jump condition (0) however specifies the post-shock pressure, so the fractional gradient is 



1 dP v 
P dij v 



4 -- (23) 



From the same jump condition we have -6 s - = ^ + 2-. For the exponential atmosphere, po — ipov/h at 



Ps Po 

the dense and tenuous ends of the remnant respectively. Thus 



Ps 


_ Po\ 






\Ps 
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or finally, at the ends of the remnant, 



ldP 
P~d^ 




(24) 



(25) 



In this equation, the leading term on the right-hand side is the positive gradient term present in the Sedov 
solution for a uniform medium. The second term is the correction factor we are looking for, positive on the 
high density end, negative on the low density end. It expresses the fact that the higher density end has a 
more rapid deceleration and therefore a steeper pressure gradient. 

For the Sedov case, equation (||) gives —Ps/Psv = 3/R. With the semimajor axis a in place of R, we 
use equation (E5h to estimate the relative gradients between the two ends at 



1 dP 



' / T>cnsc 3h 



\P dr, J 



(26) 



Sh 

Tenuous 



Upon integrating from the central pressure plateau up the slope along the major axis, the ratios of the 
post shock pressures should be similar to, but less extreme than, this ratio of gradients. As a rough 
parameterization, we can write it as 

p 1 + — 

' .S. Dense Qh (07\ 

p ~ -i a ' \ ZI > 

Tenuous n h 

where Q is somewhat greater than 3. At the end of this section, we provide a way of estimating the value 
of Q from numerical studies. 

Now, we can see that the Kompaneets approximation of constant post-shock pressure holds only for 
a < h, i.e. at times well before blowout. Our approximation of the shock by an ellipsoid is indistinguishable 
from the Kompaneets solution at this stage (the Arj/rj deviations smaller than 10~ 3 ). 

The post-shock pressure can be expressed as a function of time only, when a in equation (^7|) is 
approximated by the ID Sedov solution 

P s <x(l±(i/iQ) 2/5 )i- 6/5 , (28) 

where the top sign refers to the dense end, the bottom to the tenuous one, and (t/tg) 2 ^ 5 has replaced 
a/Qh. Keeping in mind that we are looking for a first order correction, we used the fact that at early times 
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Ps oc t 6 / 5 . For a strong, non-radiative shock, v ~ \/Ps/pa oc e z ^ 2h y / l ± {t/tq) 2 / 5 t 3 / 5 . Taking into 
account that = ^fp, but vh — — we can write 



- e - ZH/2h dz H oc y/l + (t/t Q )V 5 t~ 3 / 5 dt ; e~ ZL/2h dz L oc ^1 - (t/tq) 2 / 5 i~ 3/5 dt, (29) 
which after integrating gives 

e-«/ a& - 1 _ (1 + (t/t Q ) 2 / 5 ) 3/2 - 1 1 + Kt/tQ) 2 ' 8 1 + 



1-e 



^ /2 ' 1 l-(l-(t/t Q ) 2 /5) 3 / 2 l-i(*/<«) 2/8 



(30) 



where we first expanded the numerator and denominator to the two leading terms, and then substituted 
a from the Sedov solution. Note that for uniform post-shock pressure, the ratio (^) is equal 1, and after 
substituting zh = s — a, zl — s + a, one gets the formula (||) for the center shift s. 

Similarly, in the case of nonuniform post-shock pressure, we express Zl and zh in terms of a and s and 
after simple algebraic manipulations we get the final formula for the center shift correction 

exp(s/2h) = cosh(a/2/i) - j^sinh(a/2/i) . (31) 

The first order s correction is about 4 times smaller than the Ps correction in equation (p7|). The center 
shift correction derived above applies to models with no cooling; for numerical models, one could plot the 
quantity 

cosh(a/2/i) — exp(s/2/i) 1 



(a/2h)smh(a/2h) 2Q 
to evaluate the free parameter Q in our method. 



(32) 



5. Applications to hydrodynamical models and conclusions 

Assuming that the shock wave propagating in an exponentially stratified medium takes the shape of an 
ellipsoid, we were able to find explicit expressions for its size and expansion velocity as functions of time, 
and for the cooling time at both ends and on the equator. As presented in Table 1, these expressions take 
forms similar to the Sedov (1959) solution for a uniform ambient medium. Nevertheless, while the Sedov 
solution is self-similar, there are no such solutions for nonuniform media. The reason is that there are at 
least three independent dimensional parameters (in our case E, p* and h). A successful model of a SNR 
should adopt parameters that reproduce the observed quantities: for example the size, expansion velocity 
and progress of shell formation. If a pulsar is seen, the age can be derived as well. Our simple formulae 
provide quick ways to explore this 3-dimensional parameter space. 

A useful approach is to generate topographic (or contour) plots on the p*-h plane for assumed values 
for E and linear size of the remnant. The contoured quantities include the time required to reach the 
assumed size, the post-shock pressure, the cooling and shell formation time-scales for the dense, equatorial, 
and tenuous directions, and the shock velocities in those directions. Although the Kompancets model, like 
the Sedov model, cannot be used very reliably to estimate quantities after shell formation has occurred, we 
make the usual approximation that shortly after shell formation the shell velocity is about three-quarters 
of the shock velocity predicted for the non-radiative evolution. By this means of exploring the parameter 
space, a set of parameters fully reproducing the observables for SNR W44 was found and is given in Paper 
II. 
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There are other asymmetric SNRs which can be interpreted in terms of expansion into a nonuniform 
medium, for example a partially open HI shell in CTB 80 (Koo et al.1990), and an incomplete shell of CTA 
1 (Pineault et al.1993). The semi-circular shape of CTB 109 is ascribed to interaction with a molecular 
cloud (Tatcmatsu et al.1990). G 84.2 -0.8 (Feldt & Green 1993) shows a striking resemblance to W44 in 
radio continuum, though observations of the H I shell are not convincing. 

In the catalogue of SNRs with HI shell emission (Koo & Hciles 1991, Table 3) all sources show 
high-velocity HI gas on one side only: either receding or approaching. Heiles (1979), who catalogued 
HI shells in our Galaxy, noticed that for most of them, he could see only the approaching or only the 
receding hemisphere - a fact that he called "disturbing". Putting aside possible observational biases, a 
partially formed shell appears to be more a rule than an exception. We want to point out that explanation 
of such asymmetric SNRs or shells does not require any abrupt density change, such as encountering the 
edge of a molecular cloud (Dohm-Palmer & Jones 1996). The method of estimation of the cooling time 
presented in §3 shows that the ratio of the cooling times on the dense and tenuous ends of the remnant is 
of the order of the ratio of ambient density extremes. We can see one side of the HI shell only, because at 
first there is only one side, and then by the time the other side is formed, either the part expanding into 
denser medium has slowed down and its emission blends with that of local gas or the tenuous end remains 
invisible due to the high column density contrast between low and high density ends (Silich 1992). 

Our explicit analytical approximation provides an efficient tool for a quick exploration of the space 
of initial parameters for the SNR models. It can be used to select the initial parameters of hydro runs, 
which in turn can give us a detailed insight into the structure of the remnant, and can verify the above 
conclusions. 

We would like to thank the referee, S.A. Silich, for valuable comments, and Linda Sparke for helpful 
suggestions on the paper. DPC was supported in part by NASA Grant NAG5-3155 to the University of 
Wisconsin. 
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Table 1. 

Approximate Formulae for a SNR Expanding in a Stratified Medium 



SEDOV ID ELLIPSOIDAL 2D (EXPONENTIAL STRATIFICATION) 

' - r " 3 © ~" 2 rW2 - •"- ' = ( ~ s " (£) " /2 (2 " )S/2 1 (4) 

r 5/2 n 1/2 f2fc VV 2 n 1/2 / 

* H = 17.32 t [ yr ] = 17.32 ^ ± I (• 



E\[ 2 " Elf V2fc, 

where /(a;) = - / dy — „ arctan sinh?/ ~ =-r- for x < 3 

2i ' cosh 2 y 3 + 2a; 5 /2 

„. = ^ (^) 1/2 r^ velocity -size „. = |^ (ff ^ U (£) 

v s [km/s] = 22600 (^j ^ r pc 3 / 2 v s [km/s] = 22600 ^ (2h pc )~ 3 / 2 f ± (|-) 

where /±(l ) = -2^L22^f_ • 
Vx arctan smhx 



£ 3/14 f2/i ) 5 / 2 n 1/2 

tcooi [yr] = 4.6 x 10 4 -^j=- cooling time t cool [yr] = 17.32 1 pc \ /2 * [l(xf) + Ag±(xf)] 
nj E 5 { 



h — C ° S ^ X ex P("F^ :E ) * ^4 6 2 x 10 8 1 

W ^ ± X x 3 / 2 (arctan sinh x) 3 ' ' nl(2h pc ) 7 

, i 5 x 2 exp(±5x) (arctan sinh x) 4 

and x* is set by A = - 



2 cosh 6 x (±5eT* + ^coshx + 3 \ 

V 2a; arctan sinh x / 



Units: /i pc [pc], r pc [pc], n* [cm 3 ], E 51 [10 51 erg] 

* the top signs refer to expansion into the dense medium, the bottom into the tenuous medium; for the shock expansion velocity in the 
equatorial direction, the exponential factor should be dropped from f± (x) 



